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In this review, we show how advances in the theory of magnetic pseudodifferential op- 
erators (magnetic 'I'DO) can be put to good use in space-adiabatic perturbation theory 
(SAPT). As a particular example, we extend results of 1241 to a more general class of 
magnetic fields: we consider a single particle moving in a periodic potential which is 
subjectd to a weak and slowly-varying electromagnetic field. In addition to the semi- 
classical parameter e 1 which quantifies the separation of spatial scales, we explore 
the influence of an additional parameter A that allows us to selectively switch off the 
magnetic fleld. 

We find that even in the case of magnetic fields with components in (K''), e. g. for 
constant magnetic fields, the results of Panati, Spohn and Teufel hold, i. e. to each 
isolated family of Bloch bands, there exists an associated almost invariant subspace of 
L'^{W^) and an effective hamiltonian which generates the dynamics within this almost 
invariant subspace. In case of an isolated non-degenerate Bloch band, the full quantum 
dynamics can be approximated by the hamiltonian flow associated to the semiclassical 
equations of motion found in 1241 . 
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1. Introduction 

A fundamental and well-studied problem is that of a Bloch electron subjected to 
an electric field and a constant magnetic field where the dynamics is generated by 

H = H{e,X) -.^ ^{^iVa:~ XA{ex)y + Vr{x) +4){ex) (1.1) 

acting on L^(R!^). Here, e and A are dimensionless nonnegative parameters whose 
significance will be discussed momentarily. A is assumed to be a smooth polyno- 
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2 G. De Nittis and M. Lein 

mially bounded vector potential to the magnetic field B = dA which is constant in 
time and uniform in space. Since the magnetic field is uniform, A needs to grow at 
least linearly. The potential generated by the nuclei and all other electrons Vr is 
periodic with respect to the crystal lattice |3I4| 

r := {7 e R'' I 7 = EU ' ^ ^} (1-2) 

where the family of vectors {ei, . . . , Cd} which defines the lattice forms a basis of 
W^. The potential assumed to be infinitesimally bounded with respect to —^A^- 
By Theorem XIII. 96 in this is ensured by the following 

Assumption 1.1 (Periodic potential). We assume that Vr is F-periodic, 
i. e. Vr(- + 7) = Vr{-) for all 7 e F, and /^^ dy \Vr{y)\ < 00. 

The dual lattice F* is spanned by the dual basis {e*, . . . ,ej}, i. e. the set of 
vectors which satisfy ej ■ = 2n6kj ■ The assumption on Vr ensures the unperturbed 
periodic hamiltonian 

Hpe. = U-^y.f + Vr (1.3) 

defines a selfadjoint operator on the second Sobolev space H'^{M.'^) and gives rise to 
Bloch bands in the usual manner (cf. Section [^T|) : the unitary Blocli-Floquet-Zak 
transform Z defined by equation (j2.ip decomposes Hp^y into the fibered operator 

H^,, := ZHp,,Z-' = r dkH^,,{k) ^ dk (i(-*V, + kf + Vr{y)) 
where we have introduced the Brillouin zone 

M* {fc e M'* I fc = Yf^^i ajfij, e ["V2, +V2]} (1-4) 
as fundamental cell in reciprocal space. For each k G M* , the eigenvalue equation 

Hp,,{k)^n{k) = En{k) lfn{k), (^„(fc) £ {T^) , 

where := K'^/F, is solved by the Bloch function associated to the nth band. 
Assume for simplicity we are given a band which does not intersect or merge 
with other bands (i. e. there is a local gap in the sense of Assumption 13. 1|) . Then 
common lore is that transitions to other bands are exponentially suppressed and 
the effective dynamics for an initial state localized in the eigenspace associated to 
E^, is generated by E^,{~iVx) |6|lj . 

If we switch on a constant magnetic field, no matter how weak, the Bloch bands 
are gone as there is no BFZ decomposition with respect to F for hamiltonian p.ip . 
As a matter of fact, the spectrum of if is a Cantor set [7] if the flux through the 
Wigner-Seitz cell 

M [yeR'\y^ E -=1 "^'Sj' ^ t^V^, +V2]} (1-5) 
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Applications of Magnetic '^DO Techniques to SAPT 3 

is irrational. Even if the flux through the unit cefl is rational, we recover only 
magnetic Bloch bands that are associated to a larger lattice F' D F. A natural 
question is if it is at all possible to see signatures of nonmagnetic Bloch bands if 
the applied magnetic field is weak? 

Our main result. Theorem 14. H answers this question in the positive in the fol- 
lowing sense: if the electromagnetic field varies on the macroscopic level, i. e. e ^ 1, 
then to leading order the dynamics is still generated by the so called Peierls sub- 
stitution Eif{—iSIx — AA(ex)) + (f){ex) (defined as a magnetic pseudodifferential 
operator through equation (|2.15p , cf. Section 12. 2p . Hence, the dynamics are domi- 
nated by the Bloch bands even in the presence of a weak, but constant magnetic 
field. Furthermore, we can derive corrections to any order in e in terms of Bloch 
bands, Bloch functions, the magnetic field and the electric potential. We do not 
need to choose a "nice" vector potential for B, in fact, in all of the calculations only 
the magnetic field B enters. Existing theory is ill-equipped to deal with constant 
or even more general magnetic fields. We tackle this obstacle by incorporating the 
rich theory of magnetic Weyl calculus |20|14| with two small parameters [17] into 
space-adiabatic perturbation theory |25|24|27] . 

Magnetic Weyl calculus does not single out constant magnetic fields, in fact we 
will only assume that the components of the magnetic field B are bounded with 
bounded derivatives to any order. 

Assumption 1.2 (Electromagnetic fields). We assume that the components 
of the external (macroscopic) magnetic fields B and the electric potential </> are 
Cj^(K'^) functions, i. e. smooth, bounded functions with bounded derivatives to any 
order. 

Remark 1.1. All vector potentials A associated to magnetic fields B = dA with 
components in C{J°(M'^) are always assumed to have components in C^[(M''). This 
is always possible as one could pick the transversal gauge, 

n „i 

Ak{x):=~y2 ds Bkj{sx) sxj. 

This is to be contrasted with the original work of Panati, Spohn and Teufel where 
the vector potential had to have components in C^{M.'^). 

Under Assumptions 11.11 and II ■2[ H defines an essentially selfadjoint operator on 

The extension of the rigorous derivation of the Peierls substitution to the case 
of constant magnetic field is not only an "accademic result," it is crucial for mod- 
eling the Quantum Hall effect (QHE): the Peierls substitution can be used to link 
the QHE to the well studied Harper equation |8|2|9|11(T0] . In particular the Peierls 
substitution led Hofstadter to study a simplified tight binding model for the QHE 
today called the Hofstadter model [T^]. The Hofstadter model is a paradigm in 
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the study of fractal spectra (Hofstadter butterfly) and was used by Thouless et al. 
in the seminal paper [J^ to give the first theoretical explanation of the topologi- 
cal quantization of the QHE. More recently Avron |^ interpreted the results by 
Thouless et al. from the viewpoint of thermodynamics and connected the QHE to 
anomalous phase transition diagrams (colored quantum butterflies). Even though 
this list of publications is very much incomplete, it shows the significance of the 
Peierls substitution in the study of the QHE. The main merit of this paper is to 
provide the first rigorous proof of the Peierls substitution under the conditions rel- 
evant to the QHE and consequently the first rigorous justification for the use of the 
Hofstadter model as weak field limit for the analysis of the QHE. 

Let us now explain why we have chosen to include the additional parameter A in 
the hamiltonian H = i/(e. A). Our proposal is to model an experimental setup that 
applies an external, i. e. macroscopic electric and magnetic field. The parameter 
e ^ 1 relates the microscopic scale as given by the crystal lattice to the scale on 
which the external fields vary. We always assume e to be small. It is quite easy to 
fathom an apparaturs where electric and magnetic field can be regulated separately 
by, say, two dials. We are interested in the case where we can selectively switch 
off the magnetic field. If we regulate the strength of the magnetic field by varying 
the relative amplitude A < 1 which quantifies the ratio between scaled electric and 
magnetic field, 

B'^^ix) eXB{ex), E%x) := eE{ex). 

We emphasize that A need not be small, as a matter of fact, A = 1 is perfectly 
admissible. In this sense, A is a perturbative parameter which allows us to take 
the limit B^^'^ without changing the external electric field E^. This is very 
much relevant to experiments since magnetic fields are typically much weaker than 
electric fields and thus both, from a physics and a mathematics perspective, the 
study of the dynamics under the A — >■ limit is an interesting problem which merits 
further research. 

The aim of this review is to show how recent advances in the theory of magnetic 
pseudodifferential operators (magnetic ^'DOs) (see |14|17| as well as section [2?2| can 
be used to extend the range of validity of the results of Panati, Spohn and Teufel 
derived via space-adiabatic perturbation theory (SAPT) [24] to magnetic fields 
with components in C^{M.'^). The original proof uses standard pseudodifferential 
techniques and thus is limited to magnetic vector potentials of class C^{M.'^). In a 
recent work by one of the authors with Panati [5], adiabatic decoupling for the Bloch 
electron has been proven in the case of constant magnetic field. Their proof rests 
on a particular choice of gauge, namely the symmetric gauge. We take a different 
route: according to the philosophy of magnetic Weyl calculus, it is the properties of 
the magnetic field and not those of the vector potentials which enter the hypotheses 
of theorems. 
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Applications of Magnetic '^DO Techniques to SAPT 5 

As the proofs in |24| carry over mutatis mutandis^ we feel it is more appropriate 
to elucidate the structure of the problem and mention the necessary modifications 
in proofs when necessary. 

Our paper is divided in 5 sections: in Section [2l we will decompose the hamil- 
tonian using the Bloch-Floquet-Zak transform and rewrite it as magnetic Weyl 
quantization of an operator-valued function. Section [3] contains a comprehensive 
description of our technique of choice, SAPT. The main results, adiabatic decou- 
pling to all orders and a semiclassical limit, will be stated and proven in Section [4] 
andHJ 

Acknowledgements The authors would like to thank M. Mantoiu, G. Panati 
and H. Spohn for useful discussions. M. L. thanks G. Panati for initiating the 
scientific collaboration with M. Mantoiu. Furthermore, the authors have found the 
suggestions by one of the referees very useful. 

2. Rewriting the problem 

As a preliminary step, we will rewrite the problem: first, we extract the Bloch 
band picture via the BFZ transform and then we reinterpret the BFZ-transformed 
hamiltonian as magnetic quantization of an operator-valued symbol. We insist we 
only rephrase the problem, no additional assumptions are introduced. 

2.1. The Bloch-Floquet-Zak transform 

Usually, one would exploit lattice periodicity by going to the Fourier basis: each 
£ S{Ri) C L'^iM.i) is mapped onto 

and the corresponding representation is usually called Bloch-Floquet representation 
(see e. g. [TC]). It is easily checked that 

{J'^KKy - 7) = e-''-'{J^^){k, y) V7 G P 

holds and 7^4* can be written as 

(.F*)(fc,y) = e*^u(fc,2/) 

where u{k,y) is P-periodic in y and P*-periodic up to a phase in k. For technical 
reasons, we prefer to use a variant of the Bloch-Floquet transform introduced by 
Zak [29] which maps * e S(Ri) onto u, 

(Z*)(A:, y) := J] e-'^-(^+^) *(y + 7). (2.1) 
7er 
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The BFZ transform has the foUowmg periodicity properties: 

(Z*)(A;-7*,2/) = e+'^*-M2*)(fc,y)-:T(7*)(ZvI/)(fc,y) V7* £ T* (2.2) 
{Z^){k,y-^)^{Z^){k,y) V7Gr 

r is a unitary representation of the group of dual lattice translations F* . By density, 
Z immediately extends to i^(]RlJ) and it maps it unitarily onto 

:={^'GLL(Kfc,^'(T^)) |VXfc-7*)=T(7*)V'(fc)a. e.V7*eF*}, (2.3) 
which is equipped with the scalar product 

It is obvious from the definition that the left-hand side does not depend on the 
choice of the unit cell M* in reciprocal space. The BFZ representation of momentum 
— iVa; and position operator x on L^(MI^), equipped with the obvious domains, can 
be computed directly, 

Z{-i\7^)Z-'^ = idL2(M.) ® (-iVy) + ® idL2(T^) = -i'^y + k (2.4) 
ZxZ^^ = iVl, 

where we have used the identification Kr = L'^{M*) ® L^(T^). The superscript 
r on iVl indicates that the operator's domain Ht H H^^^^R'^ , L"^ {Ty)) consists of 
r-equivariant functions. The BFZ transformed domain for momentum ^iS7y + fc is 
L'^{M*) eg) H^{Ty). Since the phase factor r depends on y, the BFZ transform of x 
does not factor — unless we consider F-periodic functions, then we have 

ZVr{x)Z-^ = idi.(M.) Vriij) = Vr{y). 

Equations (|2.4p immediately give us the BFZ transform of H, namely 

:= ZHZ-^ = ^{-iVy + k- XAiieVl))^ + Vriy) + 0(zeV^), (2.5) 

which defines an essentially selfadjoint operator on ZC5"(M!^). If the external elec- 
tromagnetic field vanishes, the hamiltonian 

H^,, := ZHp,,Z-' = r dfc H^,,{k) (2.6) 

fibers into a family of operators on L^(T^) indexed by crystal momentum k G M* . 
T-equi variance relates Hp^j.{k — 7*) and H^^^{k) via 

H^.r{k - 1*) = r(7*) H^^^ik) r(7*)-i V7* G F* 

which, among other things, ensures that Bloch bands {En}n&i, i- e. the solutions 
to the eigenvalue equation 

H^,,{k)^n{k) = E„{k) ^„(fc), ^„(fc) e L^T'y), 
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are r*-periodic functions. Standard arguments show that H^^^{k) has purely dis- 
crete spectrum for ah k G M* and if Bloch bands are ordered by magnitude, 
they are smooth functions away from band crossings. Similarly, the Bloch functions 
k I— >■ ifnik) are smooth if the associated energy band En does not intersect with or 
touch others |26j . 

The next subsection shows that the effect of introducing an external electromag- 
netic field can be interpreted as "replacing" the direct integral with the magnetic 
quantization of H^^j. + cj). 

2.2. Magnetic ^DO and Weyl calculus 

Instead of using regular Weyl calculus, we use a more sophisticated Weyl calculus 
that is adapted to magnetic problems. It has first been proposed by Miiller in 1999 
[l9] in a non-rigorous fashion. Independently, Mantoiu and Purice [20] as well as 
Iftimie, Mantoiu and Purice [H] have laid the mathematical foundation. All of 
the main results of ordinary Weyl calculus have been transcribed to the magnetic 
context; for details, we refer to the two aforementioned publications and those we 
give in the remainder of this section. 

2.2.1. Ordinary magnetic Weyl calculus 

The basic building blocks of magnetic pseudodifferential operators are 

= P^;, := -iV, - XA{Q) (2.7) 
Q = Q, := ex. 

With this notation, H can be written in terms of P'^ and Q as the quantization of 

H{x,0-^ke + Vr{-/e) + ct>{x), (2.8) 

i. e. H = H{Q, P^). As H is the sum of a contribution quadratic in momentum and 
a contribution depending only on x, this prescription is unambiguous. However, not 
all objects (e. g. resolvents and projections) we will encounter are of this type. We 
need a functional calculus for the non-commuting family of operators Q and P"^ of 
noncomniutative observables that are characterized by the commutation relations 

z[Q,,Qj]=0 i[Qi,Pf]^e5i, t[Pf,Pf] = eXBij{Q). (2.9) 

The commutation relations can be rigorously implemented via the Weyl system 

W^{x,^) := e-'t^-Q-^-P"") =: e-^-^^^^-f^-^Q-f"*" 

where ^ is a smooth, polynomially bounded vector potential associated to a mag- 
netic field with components in C|^(R'^) and a[{x,(,), {y,ri)) := ^ • y — .t • 77 is the 
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non-magnetic symplectic form. We will also introduce the symplectic Fourier trans- 
form 

(J-./)(x,0 f dy f dr/e-"-'^)^^^^'')) f{y,r,), f e 5(r*M^), 

(271")'* jRd jRd 

which is also its own inverse on S{T*M.'^) and extends to a continuous bijection 
on the space of tempered distributions S'{T*W^). The Weyl quantization oi h € 
S{T*Ri) given by 

Op^W-tAi/ da:/ d^{T,h){x,OW^ix,0 (2.10) 

(271-)'' jRd jRd 

defines a magnetic 5*00 . Associated to this, we have a product jj^ akin to the usual 
Moyal product which emulates the product of magnetic operators on the level of 
functions on phase space, Op'^(/tJ-^5) = 0p^(/)0p^(5). For suitable functions 
/, g : r*R5^ — > C, e. g. Hormander-class symbols, their magnetic product is given 
by the oscillatory integral 

iff 9)i^^ = 77^21 f dy f dv f dz I dCe"^«-^«)'(^''')+(^'«))e'4'^«^''')'(^'';))- 
1^""^ jk;J Jr^ Jr^ jr^ 

■ e— (-^.g)(^,C) (2.11) 

where j^{x,y, z) is the scaled magnetic flux through a triangle whose corners de- 
pend on X, y and z. In |17j it was shown that for Hormander-class symbols, this 
product has an asymptotic development in e and A, 

ff9^f9-e^,{Lg}xB + Oie^), (2.12) 

where 

d d 

{f,9}\B ■■= ^{d^Jd^,g - d^Jd^,g) - A ^ By d^Jd^,9 

1=1 L] = l 

is the magnetic Poisson bracket. In the limit A — > it reduces to the standard 
Poisson bracket. The crucial fact that this product depends on the magnetic field 
B rather than the vector potential A can be traced back to the gauge- covariance 
of magnetic Weyl quantization: ii A' = A + dx is an equivalent vector potential, 
dA = B = dA' , then the quantizations with respect to either vector potential are 
unitarily equivalent, 

Op^+'^^{h) = e+^^^^f'^) Op^ih) e-'^'^('5) (2.13) 

This is generically false if we quantize /i^(a;,^) := — AA(.t)) via usual, non- 

magnetic Weyl quantization, Op^(/i) := Op(/ia).^ Otherwise, the properties of 



''Coincidentally, the quantization of polynomials of degree < 2 in momentum with respect to Op^j 
are covariant. Bloch bands, however, are not quadratic functions - as are all the other objects 
(terms of the development of the projection and the unitary) involved in this paper. 
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the relevant pseudodifFerential operators were to depend on the choice of gauge. 
Fortunately, the difference starts to appear at second order in e [T71 Section 1.1.2], 
so we generically expect results derived by nonmagnetic Weyl calculus to agree to 
first order. However, there is a second advantage of magnetic Weyl calculus (beyond 
being more natural): we can treat more general magnetic fields as properties of B 
enter rather than those oi A — and associated vector potentials are always worse 
behaved than the magnetic field. With some effort, one can treat magnetic fields 
which admit a vector potentials whose derivatives are all bounded, |9°^(a;)| < 
Ca for all a G Nq, \a\ > 1, see [TU Section 5], for instance, but with magnetic 
Weyl calculus we are instantly able to treat magnetic fields whose components are 
^oo^-j^d-j functions with zero extra effort.. The limitation to this class of fields is due 
to the fact that we are interested in Hormander class symbols which need to be 
bounded in the x variable. In fact, the extension of the results of Panati, Spohn 
and Teufel to magnetic fields of class is our main motivation as explained in 
the introduction. Covariance ensures that our results do not depend on the choice 
of a nice or symmetric gauge, any vector potential that is a smooth, polynomially 
bounded function will do. 

Many standard results of pseudodifferential theory have been transcribed to the 
magnetic case for a large class of magnetic fields: typically, it is either assumed that 
the components of B are C^i(M'*) or C{^(R'^) functions, although we shall always 
assume the latter. The quantization and dequantization have been extended to tem- 
pered distributions and it was shown that smooth, uniformly polynomially bounded 
functions on phase space Cj^i„(r*R'') are among those with good composition prop- 
erties pPl Proposition 23]. Hormander symbols are preserved under the magnetic 
Weyl product and quantizations of real- valued, elliptic Hormander symbols of pos- 
itive order m define selfadjoint operators on the mth magnetic Sobolev space |14) . 
A magnetic version of the Calderon-Vaillancourt theorem [T3] and commutator 
criteria |15| show the interplay between properties of magnetic pseudodifferential 
operators and their associated symbols. 

Lastly, we mention something that will be important in the next section: the 
magnetic Weyl quantization is a position representation for a magnetic pseudodif- 
ferential operator. However, equivalently, we can use the momentum representation 
where x is quantized to Q' := zeVj and ^ to P' := $ - \A{ieS/^). The commu- 
tation relations of the building block operators in momentum representation are 
again encoded into the Weyl system 

which is related to the Weyl system in momentum representation via the Fourier 
transform 5 : _L^(K!^) — > L^(IR^). If Op'"^ is the quantization associated to the 
Weyl system W'^, then Op'^ and Op"^ are related via ^ as well. An important 
consequence is that the Weyl product is independent of the choice of representation: 
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for suitable distributions / and g, we conclude that the Weyl products must agree: 

= d op^(/) d~':s op^(5) r'^d op^inf^^g) 

= Op'\nl,g) 

This is related to an algebraic point of view proposed in pT] and elaborated upon in 
|18| : the choice of Hilbert space in this case can be seen as a choice of (equivalent) 
representation of more fundamental C*-algebras of distributions. Properties such as 
boundedness of the quantization of certain classes of distributions and the form of 
the composition law are preserved if we choose a unitarily equivalent representation. 
In fact, we could have replaced in the above argument by any other unitary 
operator U : L^(R|^) — > % where the target space T-L is again a separable Hilbert 
space. Gauge transformations are another particular example of unitarities which 
connect equivalent representations. 

2.2.2. Equivariant magnetic Weyl calculus 

For technical reasons, we must adapt magnetic Weyl calculus to deal with equivari- 
ant, unbounded operator-valued functions. We follow the general strategy outlined 
in [21], but we need to be more careful as the roles of Q and P"* are not inter- 
changeable if B 7^ 0. We would like to reuse results for Weyl calculus on T*W^ - in 
particular, the two-parameter expansion of the product (equation p.l2p ). Consider 
the building block kinetic operators macroscopic position R and magnetic crystal 
momentum K"^, 

R = ieVfc ® idL2(T5;) = i^^k (2-14) 
= fc - AA(R), 

in momentum representation: they define selfadjoint operators whose domains are 
dense in i^, (M^, i^(T^)) where r' stands for either r : 7* e"'^* ^' or 1 : 7* (->■ 1. 
The elements of this Hilbert space can be considered as vector- valued tempered dis- 
tributions with special properties as L^, (R^, L^(T^)) can be continuously embedded 
into 5'(R^,L2(T^)). For simplicity, let us ignore questions of domains and assume 
that h e C^{T*M.i,B{L^{T'^y))) is a hounded operator-valued function. Then its 
magnetic Weyl quantization 

Dp^(/i) :- -i-j / dr / dk{F,h){k,r)W^{k,r) (2.15) 

defines a continuous operator from iS(]R^, L^(T^)) to itself which has a continuous 
extension as an operator from iS'(M^, i^(T^)) to itself [2D1 Proposition 21]. Here, 
the corresponding Weyl system 

W^{k,r) := e-^-«'-''=)'(^'^"» ® idi.(T.) ^ e-*(^-^"'-^") 
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is defined in terms of tlie building block operators K'^ and R and acts trivially 
on L^(T^). The Weyl product /jj'^g of two suitable distributions associated to the 
quantization Dp is also given by a suitable reinterpretation of equation (j2.1ip 
as / and g are now operator-valued functions. Furthermore, we can also develop 
/tt^g asymptotically in e and A [171 Theorem 1.1]. To see this, we remark that the 
difference between the products associated to Op'^ and Dp"^ is two-fold: first of 
all, Op"^ is a position representation while £)p^ is a momentum representation. Let 
£)p'"^ be the magnetic Weyl quantization defined with respect to R' := S^~^RS' = 
er and K'"* := S'^^K'^^J = — iV^ — \A{ef), i. e. the position representation. As 
explained at the end of the previous subsection, unitarily equivalent representations, 
here Dp"^ and £)p'"^, have the same Weyl product. 

Secondly, the functions which are to be quantized by Op^ and Dp^ take values 
in C and the bounded operators on L^(Ty), respectively. The interested reader 
may check the proofs regarding the various properties of the product (J^ in |20I14| 
and |17| can be generalized to accommodate operator-valued functions, including 
Hormander symbols. 

Definition 2.1 (Hormander symbols H2))). Let m G M, p € [0,1] 

and Hi, 'H2 be separable Hilbert spaces. Then a function / is said to be in 
iS™ (i3('Hi, ^2)) if and only if for aU a, a G Ng the seminorms 



a\n—ni 



are finite where || •||/3(-Hi.-h,) denotes the operator norm on B(l-Li, 112) ■ In case p — 1, 
one also writes iS"' := 5™ 

Hormander symbols which have an expansion in e that is uniform in the small 
parameter are called semiclassical. 

Definition 2.2 (Semiclassical symbols A5™(i3(Hi,H2))). A map / : 
[0,£o) — > <5™, e /e is called a semiclassical symbol of order m G M and weight 
p G [0, 1], that is / G A<S™, if there exists a sequence {/,i},igNoi fn G <S™~"'', such 
that for all iV G No, one has 

uniformly in e in the sense that for any G No and a, a G Nq, there exist constants 
CNaa > such that 



N-l 



fs ^ ^ £ fn 



n=0 



m — N p^aa 



holds for all £ G [0, £o)- If P = 1, then one abbreviates AiS™ with AS"^. 
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Lastly, we will need the notion of r-equivariant symbols. 

Definition 2.3 (r-equivariant symbols AC^ (6(^1, 7^2))) • Let : T* — > 
U{'Hj), j = 1,2, be unitary ^-representations of the group F*. Then / e AiSq is 
T-equivariant, i. e. an element of AC^ [B{Hi,H2)) , if and only if 

/(fc-7*,r)=T2(7*)/(A;,r)Ti(7*)-^ 
holds for aU keR'/.^reWf. and 7* e F*. 

Now the reader is in a position to translate the results derived in Appendix B of 
[?7] to the context of magnetic Weyl calculus, the modifications are straightforward 
and all the necessary references have been given in this section. 

3. The magnetic Bloch electron as a space-adiabatic problem 

Our tool of choice to derive effective dynamics is space-adiabatic perturbation the- 
ory [25124127] which uses pseudodifferential techniques to derive perturbation ex- 
pansions order-by-order in a systematic fashion. We adapt their results by replac- 
ing ordinary Weyl calculus with magnetic Weyl calculus. Adiabatic decoupling only 
hinges on £ ^ 1 and does not rely on A to be small. 

3.1. The three ingredients of space-adiabatic problems 

The insight of was that the slow variation of the external electromagnetic 
field (quantified by £ 1) leads to a decoupling into slow (macroscopic) and fast 
(microscopic) degrees of freedom. This is characteristic of adiabatic systems whose 
three main features are 

(i) A distinction between slow and fast degrees of freedom: the original (physical) 
Hilbert space H = L^(Rf) is decomposed unitarily into "Hsiow ® ^fast := 
L'^{M*) ® L^(T^) in which the unperturbed hamiltonian is block diagonal 
(see diagram p.ip ). 

The names slow and fast Hilbert space are due to the operators defined on 
them: on the fast Hilbert space, the two conjugate observables are — iVy and 
y acting on Hfast = L'^i'^y), the Hilbert space associated to the Wigner Seitz 
cell AI = Ty (cf. equation p.5p ): their commutator is of 0(1). The operators R 
and K'^ (cf. equation (|2.14p ) defined on Hsiow — L?{M*) are considered slow, 
because their commutator is of 0{e). Since 'Hsiow is the Hilbert space over 
the Brillouin zone M* (cf. equation p.4p ). the dynamics of the slow variables 
(R, K"^) describes the motion across unit cells in momentum representation 
whereas the dynamics of the fast variables (y, — iVj^) describe what happens 
within the Wigner Seitz cell M . 
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(ii) A small, dimensionless parameter e that quantifies the separation of spatial 
scales. In our situation, e <C 1 relates the variation of the external electro- 
magnetic field to the microscopic scale as given by the lattice constant. In 
addition, we have the parameter A. However, only the semiclassical parameter 
e is crucial for adiabatic decoupling. 

(iii) A relevant part of the spectrum, i. e. a subset of the spectrum which is sepa- 
rated from the remainder by a gap. We are interested in the dynamics asso- 
ciated to a family of Bloch bands {i?„}„gi that does not intersect or merge 
with bands from the remainder of the spectrum. 

Assumption 3.1 (Gap condition). The spectrum of /fj^^ satisfies the gap con- 
dition, namely there exists a family of Bloch bands T- — [-^-, -^+] HNq such 
that 

^inf^ dist(|J {i?„(fc)}, |J{ii;,(fc)}) =: C, > 0. 

The spectral gap ensures that transitions from and to the relevant part of the 
spectrum are exponentially suppressed. Band crossings within the relevant part of 
the spectrum are admissible, though. 

In the original publication, an additional assumption was made on the existence 
of a smooth, r-equivariant basis, a condition that is equivalent to the triviality of 
a certain U{N) bundle over the torus where TV := \I\ is the number of bands 
including multiplicity. At least for the physically relevant cases, i. e. d < 3, Panati 
has shown that this is always possible for nonmagnetic Bloch bands [23]. For d > A, 
our results still hold if we add 

Assumption 3.2 (Smooth frame {d > 4)). If d > 4, we assume there exists 
an orthonormal basis (called smooth frame) {^j{')}j=i \x\ of whose elements are 
smooth and r-equivariant with respect to k, i. e. Lpj{- — 7*) = T{'y*)ipj(-) for all 
7* e r* and for ah j G {1, . . . , 

3.2. Rewriting the unperturbed problem: an adiabatic point of view 

Let us consider the unperturbed case, i. e. in the absence of an external electromag- 
netic field. Then the dynamics on Tir is generated by H^^^ = dk Hp^^.{k). Each 
fiber hamiltonian H^^j.{k) is an operator on the fast Hilbert space Hfast = i^(Ty). 
Then ttq = j'^^, dfc7ro(fc) is the projection onto the relevant part of the spectrum, 
where 

'^o(fc) := ^ |v3„(fc))(^„(fc)|. 

Even though the </J„(fc) may not be continuous at eigenvalue crossings, the projec- 
tion k M> 7ro(fc) is due to the spectral gap. Associated to the relevant band is a (non- 
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unique) unitary uq = J^j, dkuo{k) which "straightens" Ht into L'^{AI^) ^L^{Ty): 
for each k G M* , we define 

where Xn G L^(T^), n E I, are /iied vectors independent of k and (also non- 

unique) acts on the complement of ran7ro(A:) and is such that uq is a proper unitary. 
Even though this means uq is not unique, the specific choices of the {xn}nei and 
Uq will not enter the derivation. Then we can put all parts of the puzzle into a 
diagram: 



n 



Z-^ttqZ I 



I 

Ar 



L2(M*)(g)L2(Td) 



(3.1) 



The reference projection Hi-cf = idL2(j\/*) ® TTiof acts trivially on the first factor, 
L'^{AI*), and projects via 

N 

TTrof = \Xj){Xj \ = uo{k) no{k) u^ik) (3.2) 

onto an iV-dimensional subspace of L^(Xy). We will identify TT,-c[L'^{Ty) with 
when convenient and in this sense, we identify the range of rtj-of with the reference 
space TLr-cf := L^{M*) ® C^. 

The dynamics in the lower-right corner is generated by the effective hamiltonian 

hcS := Hrof Uo ^ifcr U*Q Href 

which reduces to En{k) if the relevant part of the spectrum consists of an isolated 
Bloch band. 



3.3. Adiabatic decoupling in the presence of external fields 

Now the question is whether a similar diagram exists even if the perturbation is 
present, i. e. if there exist a tilted projection 11, an intertwining unitary U and 
an effective hamiltonian hcs that take the place of ttq, uq and /icffo? This has 
been answered in the positive for magnetic fields that admit C{f (K'^,M'^) vector 
potentials in [53] where these objects are explicitly constructed by recursion. We 
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replace standard Weyl calculus used in the original publication with its magnetic 
variant (cf. Section 12. 2p which naturally allows for the treatment of more general 
magnetic fields with components in C^. 

The construction of 11 and U detailed in the next section is a "defect construc- 
tion" where recursion relations derived from 

= n [n,iJ^] = 

[/*[/ = id-H,, UU* ^ idL2^^(A/')(»L2„(A/) UUU* = U^ef 

relate the nth term to all previous terms. These four conditions merely characterize 
that n and U are still a projection and a unitary (first column) and adapted to the 
problem (second column). These equations can be translated via magnetic Weyl 
calculus to 

Trf ^ = 7r + 0(e°°) [tt, H^]^^ = 0{e°°) (3.3) 

ufu* = 1 + 0{e°°) = u*fu ufnfu* = ir^ff + 0{e°^) 

where 

i/^(fc,r) :=i(-jV, + fc)' + yr(y) + 0(r) (3.4) 

is the operator-valued symbol associated to the magnetic pseudodifferential opera- 
tor defined by equation ()2.5p . Note that the magnetic vector potential A does 
not enter the definition of the symbol and we do not need to impose conditions 
on A to ensure the symbol is well-behaved after minimal substitution. This 
is why magnetic Weyl calculus can be used to treat much more general magnetic 
fields. 

For technical reasons, Dp^(u) and C/, for instance, agree only up to an error 
that is arbitrarily small in e with respect to the operator norm, 

C/ = Dp^(u) + C'||.||(e°°). 

The tilted projection and intertwining unitary are now used to define the effective 
hamiltonian as the magnetic quantization of 

Kff := TTrcf ufH^fu* TTrcf 

which generates effective dynamics, i. e. for initial states in 11?^,- we can approximate 
the full time evolution in terms of e"''?^'' C*"*'*). In turn, the effective quantum 
evolution can be approximated by semiclassical dynamics. Theorem 15. 2| the main 
result of the next section, will make this statement precise. 

4. Derivation of effective quantum dynamics 

The aforementioned "defect construction" yields the tilted projection tt and the 
intertwining unitary u as asymptotic expansion in e and A. It is important that 
the decoupling is solely due to the separation of spatial scales quantified by e and 
independent of A which regulates the strength of the magnetic field. 
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4.1. The dynamics in the almost invariant subspace 

We will quickly explain how 11 and U are computed order- by-order in s and A. We 
adapt the general recipe explained in |27| to incorporate two parameters: since the 
decoupling is due to the separation of spatial scales quantified by e <C 1, we will 
order corrections in powers of e first. Expanding the magnetic Weyl product to 
zeroth order, we can check 

n„fn„=n„ + 0{e) [tt^, H^]^, = 0{e) 

uofuo* = 1 + 0{e) = un*fuo uofnofu* = ^eff + 0{e). 

Here, [770,^?"^]^^ := ttq^^H'^ — H'^^^ttq denotes the magnetic Weyl commutator. 
The asymptotic expansion of the product is key to deriving corrections in a system- 
atic manner: the 0{e) terms can be used to infer tti and u^, the subprincipal sym- 
bols. Then, one proceeds by recursion: if tt'"^ := X^ILo ^^'^ ^'"'^ X^ILo 
satisfy equations (j3.3p up to errors of order e"+^, then we can compute 7r„+i and 
Un+i- The construction of tt and u follows exactly from Lemma 3.8 and Lemma 
3.15 of |27j : it is purely algeraic and only uses that we have a recipe to expand the 
Moyal product in terms of the semiclassical parameter e. Let us define 

T,in)f^{n) _ ^(n) ^. £«+1g„^^ + C>(£"+2) (4.1) 

as projection and commutation defects and 

(4.2) 

(^(") + e"+ia„+i«o)f 7r(«+i)f (u(«) + £"+ia„+iUo)* =: + 0(e"+') 

as unitarity and intertwining defects. The diagonal part of the projection n'^^i can 
be computed from Gn+i via 

T^n+l •= -TToGn + lTTo + (1 - 7ro)G„+l(l - TTq). (4.3) 

The term 

a-n+i = -5^1+1 (4.4) 

stems from the ansatz — {an+i +6„+i)mo where a„_|_i and &„+i are symmetric 

and antisymmetric, respectively. One can solve the second equation for 

K+i = [7rrcf,S„+i] (4.5) 

where TTj-cf is the reference projection on i^(T^) given by equation p.2p . This equa- 
tion fixes only the off-diagonal part of hn+i as ■n-cciBn+iT^xci = = (1 — TTrof )i?n+i(l — 
TTrof) and in principle one is free to choose the diagonal part of This means. 
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there is a freedom that allows arbitrary unitary transformations within 7riofi^(T^) 
as well as its orthogonal complement. In general, it is not possible to solve 

[^^,^°°] =-^n+l (4.6) 

explicitly since Bloch functions at band crossings within the relevant part of the 
spectrum (which are admissible) are no longer differentiable. In any case, tt can be 
constructed locally around (fco, rg) by asymptotically expanding the Moyal resolvent 
— z)^"^^^, i. e. the symbol defined through the relations 

{H^ - z)f{H^ - z)(-i)« =1 = {H^ ^ zf-^^^f{H^ - z), 

and setting 

T,(k,r) = ^[ dz(iJ^-z)(-i)^(/c,r) + 0(e°°) (4.7) 

Jc{ko,ro) 

in a neighborhood of {ko,ro). A recent result by Iftime, Mantoiu and Purice [T5] 
suggests that under these circumstances {H^ is elliptic and selfadjoint operator- 
valued) (i?"^ — z)^^^-'^ always exists and is a Hormander symbol even in the presence 
of a magnetic field. We reckon their result extends to the case of operator-valued 
symbols, but seeing how tedious the proof is, we simply stick to the procedure used 
by Panati, Spohn and Teufel |27[ Lemma 5.17]. This construction uniquely fixes 
the tilted Moyal projection tt, but not the Moyal unitary u. 
As the two-parameter expansion of the product 

oo n 
n=0 k=Q 

contributes only finitely many terms in A for fixed power of n of e [17| , we can order 
the terms of the expansion of tt and u in powers of A, e. g. 

n 
k=0 

The magnetic Weyl product as well as its asymptotic expansion are defined in terms 
of oscillatory integrals, i. e. integrals which exist in the distributional sense. If we 
take the limit A — >■ of fi^g, we can interchange oscillatory integration and limit 
procedure jl31 p. 90] and conclude limA-s.o /tJ^5 = fig where is the usual Moyal 
product. Similarly, we can apply this reasoning to the asymptotic expansion: for 
any fixed iV e Nq, we may write the product as 

N . n s 

ffa = E ^" E (/f <?)(«..) + <+i (/, 5) 

and taking the limit A — > means only the nonmagnetic terms (/tl^S')(n,o) sur- 
vive. The remainder also behaves nicely when taking the limit as it is also just 



January 12, 2011 3:40 WSPC/INSTRUCTION FILE rmp_mag_sapt 



18 G. De Nittis and M. Lein 

another oscillatory integral and limA->o ^N+iif^ 9) exactly the remainder of the 
nonmagnetic Weyl product. 

Hence, we can now prove the main result of this paper: 

Theorem 4.1 (Effective quantum dynamics). Let Assumptions \1.2\ 
and \ 3.1\ be satisfied. Furthermore, if d > 4, we add Assumvtion \ 3/A Then there 
exist 

(i) an orthogonal projection 11 G B{Hr), 

(a) a unitary map U which intertwines T-Lt and L'^{M*) Cg) L^(T^), and 
(iii) a selfadjomt operator Op^^ih^s) £ B{L'^{M*) (g) C'^) , N := \I\ 

such that 

\\[H^ ,U]\\ = 0{e°°) (4.8) 

and 

II (e—^^ C/*e-'^°P^('''--«-)i7)n||g(^^j = 0{e°°{l + \s\)). (4.9) 
The effective hamiltonian is the magnetic quantization of the T* -periodic symbol 

oo 

/leff := 7r,ef ufH^fu* ^rcf 51 ^" ^ ^'^-^1 (^(C^^)) (4.10) 

n=0 

whose asymptotic expansion can be computed to any order in e and A. To each order 
in e, only finitely many terms in A contribute, ft-cffn = X]fc=o ^ ^cff(n,fc)- 

One deduces from equation (|4.8p and a Duhaniel argument that the unitary 
time evolution generated by iJ"^ and 11 almost commute even for macroscopic times 
t = ss and hence, up to an error of arbitrarily large order in e. the space Hi.cf'Hr is 
left invariant by the dynamics, 

(l-n)e-'^^^n = 0||.||(e°°|t|). 

Thus we call EfHr the almost invariant subspace and in a sense, it is the tilted 
"eigenspace" associated to the family of relevant bands. 

The proof of the above theorem amounts to showing (i)-(iii) separately. 

Proposition 4.1 (Tilted projection). Under the assumptions of Theorem \4.1\ 
there exists and orthogonal projection 11 G B{'Hr) such that 

[i7^,n] =0||.||(e-) (4.11) 

and n = Dp"^(7r) + ©[[.[[(e"") where £)p'^(7r) is the magnetic Weyl quantization of 
a T-equivariant semiclassical symbol 

oo 

^X^e"^„ G AC^(S(?^fast)) 
n=0 
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whose principal part TTQ{k,r) coincides with the spectral projection of {k,r) onto 
the subspace corresponding to the given isolated family of Bloch bands 
Each term in the expansion can be written as a finite sum 

n 
fc=0 

ordered by powers of X. For A — !■ 0, the projection n reduces to the nonmagnetic 
projection tt" x J2n=o ^(".o) • 

Proof. [Sketch] The proof rehes on a weU-developed magnetic Weyl calculus 
adapted to operator-valued symbols (cf. Section 12. 2p and the gap condition. In 
particular, one needs a magnetic Calderon-Vaillancourt theorem, composition and 
quantization of Hormander symbols |14) and finally, an asymptotic two-parameter 
expansion of the magnetic Weyl product jl^ [17]. The interested reader may check 
line-by-line that the original proof |27[ Proposition 5.16] can be transliterated to 
the magnetic context with obvious modifications. If we were using standard Weyl 
calculus, the major obstacle would be to control derivatives of tt since vector po- 
tentials may be unbounded. In magnetic Weyl calculus the vector potential at no 
point enters the calculuations and the assumptions on the magnetic field assure 
that TT G AC^ [B{'Hiast)) is a proper r-equivariant semiclassical Hormander-class 
symbol (cf. Definition 12. 3p . 

The fact that we can write all of the 7r„ as finite sum of terms ordered by powers 
of A stems from the fact that calculating 7r„ involves the expansion of the product 
up to nth power in e, e. g. for the projection defect, we find 

a-\-b-\-c—n 

= E E E E (-(a.aof + Oie-^')- 

a+b+c=n a'=0 b'=0 c'=0 

Certainly, the exponent of A is always bounded by n > a' + b' + c'. And since the sum 
is finite, this clearly defines a semiclassical symbol in e is shown (cf. Definition [22]) • 
Similar arguments for the commutation defect in conjunction with the comments in 
the beginning of this section show tt to be a semiclassical symbol. It is well-behaved 
under the A — > limit and reduces to the projection associated to the case B = 0. 

Lastly, to make the almost projection Dp'^(7r) into a true projection, we define 11 
to be the spectral projection onto the spectrum in the vicinity of 1, 

n := / dz(Dp^(7r) 
This concludes the proof. □ 
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Similarly, one can modify the proof of Proposition 5.18 found in |27) to show 
the existence of the intertwining unitary. 

Proposition 4.2 (Intertwining unitary). Let {En}nei be a family of bands 
separated by a gap from the others and let Assumption \1.1\ be satisfied. If d > 
3, assume Uq £ ^Sq (;B(7^fast)) • Then there exists a unitary operator U : T-Lt — > 
L^{M*)®L'^{T^) such that U ^ Dp^{u) + Oi\\{e'=^) where 

OO 

ux^£"7/„ eA5° (Belfast)) 

is right-T-covariant at any order and has principal symbol Uq. Each term in the 
expansion can be written as a finite sum 

n 
fc=0 

ordered by powers of X. For A ^ 0, the unitary u reduces to the nonmagnetic unitary 

" ^ T,n=0^ "(",0)- 

Proof. [Sketch] Equations ()4.4p and ()4.5p give us a„+i and bn+i which combine to 
Un+i = (a„+i + 6„+i)iio; by Theorem 1.1 from |17| it is also in the correct symbol 
class, namely iSq (B(Hfast))- The right r-covariance is also obvious from the ansatz. 
Lastly, the true unitary U is obtained via the Nagy formula as described in |27I|II 

Proof. [Theorem 14.11 (Sketch)] The existence of IT and U have been the subject 
of Propositions 14.11 and 14.21 By right- r-covariance of u, h^g is a r*-periodic sym- 
bol and since it is the magnetic Weyl product of {T*E.'^ , B {L'^ {T^))) functions. 
Theorem 1.1 from [T7] ensures that the product and its asymptotic two-parameter 
expansion are in C^{T*R'^,B{L'^{T^))) as weh. Equation gH) follows as usual 
from a Duhaniel argument. □ 

The crucial statement of Theorem 14.11 is equation ()4.9|) and it is worthwhile 
to discuss its implications to applications: in practice, one only computes finitely 
many terms of the asymptotic expansions oi H, U and ft-off • Let us call their finite 

resummations w(') := X)L=o '^'''^ •= J2n=o ^" '^'i ^^'^ '^cff''' := Z]L=o ^offn- 

Assume we are interested in times t = e^s where \t\ < r. Then a closer inspection 
of the Duhamel argument in the proof of Theorem 14.11 yields 

Hence, if we would like to consider macroscopic times t = es and make an error 
in the propagation of order ©(e^), then we need to expand the tilted projection, 
intertwining unitary and the effective hamiltonian to second order. However, if the 
relevant band consists of a single band, computing the first-order correction to the 
effective hamiltonian suffices as we shall see in Section [5l 
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4.2. Effective dynamics for a single band: the Peierls substitution 

In case the relevant part of the spectrum consists of a single non-degenerate band 
, we can calculate the first-order correction to explicitly: the magnetic Weyl 
product reduces to the pointwise product to zeroth order in e. Thus, we can directly 
compute 

^offo = '"'rcf Mo Hq Uq TTrcf =: TTicf ho TTrof = + (j). 

For the first order, we use the recursion formula [27l eq. (3.35)] and the fact that 
/loffo is a scalar-valued symbol: 

/loffi = (ui Hq - houi + (uott"^-ffo)(i) - (/iotJ"^wo)(i))wS 

= TTrcf [uiUQ,ho] TTrct + (uotJ"^^^o)(l)Uo - (ho'i^ Uq) {1)Uq 

The term with the magnetic Poisson bracket can be easily computed: 

TTrof (^{'^0,Ho}^g - {ho^Uaj^g^ Wq TTrcf = 

= TTrcf (dki Uo Hq - XBij (^fej UO ^ dkj Uq Uq ho - /iq Uq Ok, U^) + 

+ dri hQ dk, Uq + XBij dk, hQ dk, U(^ Uq TTrcf 
= 2i(dri4) - XBij dkjE.^.) Al + XBij TTrcf dkiUQ Uq {dkjUQ Uq ho + ho Uq OkjUg) TTrcf 

= 2z(a,> - XBij dk^E,) Al + XBij {dk,^k, {Hp,, - E,)dk^^b) 

The first term combines to a Lorentz force term, the second one - which is purely 
imaginary - yields the Rammal- Wilkinson term. The components of the magnetic 
field are Cf^(R'') functions and hence, principal and subprincipal symbol are in 
C^{T*R'^) as well. 

This means, we have proven the following corollary to Theorem 14. II 

Corollary 4.1. Under the assumptions of Theorem \4-.l\ the principal and subprin- 
cipal symbol of the effective hamiltonian for a single non-degenerate Block band E^ 
are given by 

h,SQ=E,+(j) (4.12) 
Kei^ -{-dr,(j) + XBij dk^E^) Ai - XBi^ Mij 
=: -Fi^oriAi-XBijMij. 

where 

Ai{k) ■.= i{Lph{k),Vkfh{k)) 

and 

Mij{k) = Re (§(afc,(^b, (i?pcr(fc) - E,{k)) dk^^h)) 
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are the Berry connection and the so-called Rammal- Wilkinson term, respectively. 

To leading order, this is the well-known Peierls substitution. In particular, for 
zero electric field, = 0, and constant magnetic field (and the equations written in 
symmetric gauge), the ansatz 

d 

E^{k) = cos kj 
yields the celebrated Hofstadter model p2] . 



5. Derivation of semiclassical equations of motion 

In the preceding section, we have approximated the full quantum evolution by a 
simpler effective quantum evolution on a smaller reference space Hrd = L'^{M*) 



C . Now, we will link the time evolution generated by hcs - and thus 
to a semiclassical flow which contains quantum corrections. Conceptually, this is a 
two-step process: if we reconsider the diagram of spaces. 



i2 



I 
I 

I 



n 



-^L'^{M*)®L^{Ti,) 



n,of 



(5.1) 



we notice that our physical observables live on the upper- left space L^(R^) - or 
equivalently on T-Lt-. The effective evolution generated by /loff approximates the dy- 
namics if the initial states are localized in the almost invariant subspace associated 
to the relevant bands. In this section, we always assume the relevant part of the spec- 
trum consists of a single non- degenerate hand E^, and thus L^{M*)^C^ = L'^{M*). 

In a first step, we need to connect the semiclassical dynamics in the left column 
of the diagram with those in the lower-right corner. The second, much simpler step 
is to establish an Egorov-type theorem on the reference space Hici = L^(M*) (g)C^. 



5.1. Relation between dynamics for macroscopic and effective 
observables 

Since we are concerned with the semiclassical dynamics of a particle in an electro- 
magnetic field, the magnetic field must enter in the classical equations of motion. 
There are two ways: either one uses minimal coupling, i. e. one writes down the 
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equations of motion for position r and kinetic momentum k ~ k' — XA{r) with 
respect to the usual symplectic form for the variables {k',r) and the hamiltonian 
H'^ {k' — \A{r),r). Then the classical flow which enters the Egorov theorem is 
generated by 

where the appearance of B in the matrix representation of the symplectic form 
is due to the fact that k is kinetic momentum. What constitutes a suitable ob- 
servable? Physically, we are interested in measurements on the macroscopic scale, 
i. e. the observable should be independent of the microscopic degrees of freedom. 
On the level of symbols, this means /(/c, r) has to commute pointwise with the 
hamiltonian H^{k,r) for all k and r. Hence, such an observable is a constant of 
motion with respect to the fast dynamics. In the simplest case, the observables are 
scalar- valued. This also ensures we are able to "separate" the contributions to the 
full dynamics band- by-band. Note that this by no means implies D\)^{f) commutes 
with Op^{H^), but rather that all of the non-commutativity is contained in the 
slow variables {k,r). 

Definition 5.1 (Macroscopic semiclassical observable). A macropscopic ob- 
servable / is a scalar- valued semiclassical symbol (cf. Definition [22]) AiSq (C) which 
is r*-periodic in fc, f{k + 7*, r) = f{k, r) for aU (fc, r) e T*R'^, 7* £ T*. 

Our assumption that our dynamics lives on the almost-invariant subspace IIT^t- 
modifies the classical dynamics to first order in e as well: instead of using K'^ 
and R as building block observables, the proper observables should be IIK'^II and 
nRn. Equivalently, we can switch to the reference space representation and use the 
magnetic quantization of 

k^s := TTrcf ufkfu* TTi-ct = fc + eXB{r)A{k) + 0{e^) (5.3) 
rcff := TTrcf ufrfu* 7r„f = r + eA{k) + 0{e') 

The crucial proposition we will prove next says that for suitable observables /, the 
effect of going to the effective representation is, up to errors of order at least, 
equivalent to replacing the arguments k and r by fccff and r^s, 

nreft/Dp^(/)f/*n,ef - Df^ {n,,f uf ff U* TT^cf) + 0||.||(e°°) 

=:Dp^(/eff) + 0||.||(e°°). (5.4) 

Then it follows that the effective observable /eg coincides with the original observ- 
able / after a change of variables up to errors of order 0{e'^), 

/off = 7r„f ufffu* ^,ef = / o Teff + 0{e^), (5.5) 
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where the map Toff '■ {k,r) t— >■ (fcoff,^cff) maps the observables k and r onto the 
effective observables fcofr and rofr defined via equations (|5.3p . 

Proposition 5.1. Let f be a macroscopic semiclassical observable. Then up to 
errors of order equation (|5.5p holds and consequently, we have 

n,efC/Dp^(/)C/*nrcf =Dp^(/cfi-) +0||.||(e°°) =£>P^(/°reff) +0||.||(£2). (5.6) 

Proof. The equivalence of the left-hand sides of equations (|5.6p and (|5.5p follows 
from U = Dp'^iu) + 0\^.\^{e°°) and the fact that we only need to consider the first 
two terms in the s expansion. With the help of Theorem 1.1 from |17| . we conclude 
/cff S AiS° is also a semiclassical symbol of order 0. The left-hand side of ()5.5p can 
be computed explicitly: to zeroth order, nothing changes as / commutes pointwise 
with u and u*, 

/cffO = Tl'icf Mo/Wq ""rcf = /o- 

To first order, we have 

/cffl = TTrcf (uofl +Ul/o - /cffOUl + (wotJ'^/o)(l) - (/off ott^ "o) (1) ) TTrof 

= TTrof UoflUQTTrcf " f ({"O, fo}xB ^ {fcffOi "o};^^) 
= fl~i {drj /o + ASjj (r) dk, /o) TTrcf dkj Uq Uq TTi-ef 

= .fi + idrJo + XBi,{r)dkJo)Aj. 

On the other hand, if we Taylor expand / o Teg = /(fcoff , ''off) to first order in e, we 
get 

/(fcoff, ''off) =fo{k + eXB{r)A{k) + 0{e^),r + eA{k) + 0{e^)) + 

+ efi{k + eXB{r)A{k) + 0{e^),r + eA{k) + 0{e^)) + 0{e^) 

= /o(fc,r) +e(/i(fc,r) + Aafe,/o(fc,r) B,,(r)^-(fc)+ 
+drJo{k,r)Aj{k))+0{e^) 
which coincides with /cff up to 0{e^). □ 

Now if the equations of motion (|5.2p are an approximation of the full quantum 
dynamics, what are the equations of motion with respect to the effective variables? 
The classical flow generated by h^s with respect to the magnetic symplectic 
form (equation (|5.2p ) can be rewritten in terms of effective variables, 

^macro ^ j.^^ ^ ^off ^ _^ ^(^2)^ (5 7) 

The right-hand side does not serve as a definition for the flow of the macroscopic 
observables, but it is a consequence: $™acro ^j^g g^.^ associated to a modified 
symplectic form and a modified hamiltonian. The modified symplectic form includes 
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the Berry curvature associated to acting as a pseudo-magnetic field on the 
position variables. 

Proposition 5.2. Let <i)™acro y^^^ T*W^ generated by 

hscikcS,res) (5.8) 



AB(reff ) -id W ''cff \ / V^^„. 
^ +id sn{k^s) Ucffj IVfc^,, 



where the semiclassical hamiltonian is given by 

hsc ■■= hcs o T~g . (5.9) 

Then equation (|5.7p holds, i. e. and Tcff ° o r^~ff^ a^ree up to errors of 

order . 

Proof. We express k and r in terms of fcoff and Toff in (|5.2p since, for e small enough. 
Toff : (A:,r) t-> (fcoff,roff) is a bijection. For instance, the semiclassical hamiltonian 
/icff o T~Q simplifies to 

^sc(fcGff,fGff) :=^cff(fcofr - £AS(roff)^(fcofT),roff - eA{Kfi)) + 0[e'^) 
= {E,{k,s) + 0(refr)) - eAB(reff) • M{k,s) + Oie"). 

The symplectic form can be easily expanded to 

/ AB(reff - £^(fccff )) -id\ ^ 
+id j 

^ (^B{r,s) -id^ _^ |^Aa,„„,S(refr)^(/ceff) 0^ _^ ^(^2)^ 

The other two terms, the time derivatives and gradients of fcoff and refi have slightly 
more complicated expansions, but they can be worked out explicitly. Then if we put 
all of them together, we arrive at the modified symplectic form (|5.8p . This proves 
the first claim. 

Hence, the hamiltonian vector fields agree up to 0(e2) and Lemma 5.24 in [571 
implies that also the flows differ only by 0(£^). □ 

Remark 5.1. These equations of motion have first been proposed in the appendix 
of [13] and we have derived them in a more systematic fashion. The effective coor- 
dinates Toff and /coff are associated to a non-standard Poisson structure on r*R^: 
from equation ()5.8p . one can read off that the Poisson bracket with respect to reff 
and fccff is given by 

d d 

1=1 
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and thus different components of position r^ff no longer commute, 

Hence, Q acts as a pseudomagnetic field that is due to quantum effects. 
5.2. An Egorov-type theorem 

The semiclassical approximation hinges on an Egorov-type theorem which we first 
prove on the level of effective dynamics: 

Theorem 5.1. Let hcs be the effective hamiltonian as given by Theorem \4-l\ as- 
sociated to an isolated, non-degenerate Block band E^,. Then for any T* -periodic 
semiclassical observable f G AiS", / = /o +e/i, the fl.ow generated by h^s with 
respect to the magnetic symplectic form ( equation (j5.2p ) approximates the quantum 
evolution uniformly for all t G [— T, +T], 

|e+*^°P''(''"«)Dp^(/)e-*^°P''(''"«) -Dp^(/o$f )||^^^^^^^^^^ < Ce^. (5.10) 

Proof. Since h^s G C^{T*W^), the flow inherits the smoothness and /o$^ff , A (/o 
e A5°(C) remain also r*-periodic in the momentum variable. To compare the 
two time-evlutions, we use the usual Duhammel trick which yields 

+jf Op'*(/ieff) £)p'^(/) e-*f '^P'^C'off) - Op^(/ o $f ) = 

= As e+^f ^P'C^o..) . [Dp^(/.eff ), £)p^(/ o <I>f ,)] + 

-Dp^(£(/o<,)))e-f^P-'(''-) 
= ^* As e+^f Dp^ ( J [/.eff , / o + 

-{K^,fo^f_^]\e-^i^^'(^-^\ 



e 



(5.11) 

The magnetic Moyal commutator - to first order - agrees with the magnetic Poisson 
bracket, 

I [Ks, f o $f ,]^, = / o cf>f + Oie') 

d d 

= X! (^'^i ^''^ ~ ^t^ff ^'^i /) ^ X! -^'j ^'^i ^<=ff ^'^j / + ^('^^) ■ 

Hence, the term to be quantized in equation (j5.1ip vanishes up to first order in e, 
r.h.s. of (EH]) = rdse+*t'^P''(''-«)Dp^(0 + O(£2))e-^f°P''(''<=«-) =0||.||(£2). 

"'0 
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This finishes the proof. □ 

The main result combines Proposition 15.11 with the Egorov theorem we have 
just proven: 

Theorem 5.2 (Semiclassical limit). Let Assumptions M . 1\ \3.1\ he satisfied; 
if d > 4:, assume in addition that Assumption holds true. Furthermore, let 
us assume the relevant part of the spectrum consists of a single non- degenerate 
Bloch band E^. Then for all macroscopic semiclassical observables f (Definition \5.1\) 
the full quantum evolution can be approximated by the hamiltonian flow cl)™a<=''° as 
given in Proposition \5.2\ if the initial state is localized in the corresponding almost 
invariant subspace Z~^IIZL'^{'E.'^), 



(5.12) 

Proof. We now combine all of these results to approximate the dynamics: let / 
be a macroscopic observable. Then if we start with a state in the range of Z~^nZ, 
the time-evolved observable can be written as 

Z-^nZe^'^^' Z-iDp^(/)Ze+'i^'z-inZ = 

= Z-'Ue-^i"" Dp^{f)e+'i^"uZ 

= Z-^ U-^U,,tUU-^e-'i''UO\)^{f)U-^e+'^^UU-'U,,fU Z + Oii.\iie°°) 
= Z-i U-^n,^{e-'i''UOp'^if)U-^e+'i''U,^fU Z + O,,.,, {e°°) 



Z-i C/-in,efe-*i''<=«nrcf C/Dp^(/)[/-in,.efe+'i'^"«n,.ef C/Z + 0||. II 



After replacing U with Dp^{u) (which adds another 0||.||(£°°) error) and Href 



with Dp"^(7ri.of ), the term in the middle combines to the quantization of fcs = 
""■rof wtt'^/tt^w* '"'rcf • We apply Proposition 15.11 and the Egorov theorem involving 
/leff and obtain 

. . . = Z-i U-^n,,fe-'if'^« Dp^(7r,et ufffu* Tr,,f) e+^^'^-^n^eff/ Z + 0||.||(e°°) 
= Z-i [/-in,efe-**^"« Op^if^s) e+*i^^«n,etC/Z + 0||.||(e°°) 
= Z-in[/~iDp^(/oTeff o$f ) [/nZ + C'||.||(e2). 

Since two flows are 0{£^) close if the corresponding hamiltonian vector fields are 
P71 Lemma 5.24], we conclude 

. . . = Z-i n [/-I Dp'^(/ o Teff o $f o T^-i o Teff) C/nZ + 0||.||(e2) 

= z-i n Op-^{f o $5"'^"° o Tcff ) [/ n z + Oii.ii (e°°) 
^ z-^ nDp'^(/ o $j°a"°) nz + e»||.||(£2). 

This finishes the proof. □ 
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